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, — ,! ABSTRACT 

a,; 

^ e— | Focusing of electromagnetic plane wave from a large size paraboloidal reflector, 

& : 

composed of layers of chiral and/or chiral nihility metamaterials, has been studied us- 

■ ing Maslov's method. First, the transmission and reflection of electromagnetic plane 

> '■ 

wave from two parallel layers of chiral and/or chiral nihility metamaterials are cal- 

i> : 

culated using transfer matrix method. The effects of change of angle of incidence, 

<n : 

^ ! chirality parameters and impedances of layers are noted and discussed. Special cases 

by taking very large and small values of permittivity of second layer, while assuming 
value of corresponding chirality equal to zero, are also treated. These special cases 
are equivalent to reflection from a perfect electric conductor backed chiral layer and 
nihility backed chiral layer, respectively. Results of reflection from parallel layers have 
been utilized to study focusing from a large size paraboloidal reflector. The present 
study, on focusing from a paraboloidal reflector, not only unifies several published 
works conducted by different researchers but also provides better understanding of 
new cases. 



1.0 Introduction 

Chiral media, have been well known for a long time due to some interesting 
properties, exhibiting in the optical range of frequencies. Out of these properties, 
optical activity and circular dichroism are of practical interest [1-6] . Chiral medium can 
be thought as composed of numerous randomly oriented chiral objects which can never 
be brought into congruence with their mirror images by any translation or rotation [7- 
10]. Many chiral objects are found in nature such as irregular tetrahedrons, sugar 
molecules and wire helices. In fact the word "chirality" is derived from the greek 
word "cheir" having meaning of hand; an object possessing above mentioned property. 
In science particular word "chirality" instead of dissymmetry was first introduced by 
Lord Kelvin, Professor of natural philosophy in the University of Glasgow. Now this 
word is used to describe the media microscopically composed of the chiral objects. 

When linearly polarized wave falls on a slab of chiral medium, it splits into two 
circularly polarized waves: one left circularly polarized and the other right circularly 
polarized [8-10]. After passing through the slab the two waves combine to yield a 
linearly polarized wave whose plane of polarization is rotated with respect to the plane 
of polarization of the incident wave. Effect of chirality of the medium which rotates 
the plane of linearly polarized wave passing through it is studied by Biot, Arago 
and Friesnel. This phenomenon is named as optical activity. The amount of rotation 
depends upon the distance traveled by the wave in the slab. This phenomenon is named 
as circular birefringence. Circular dichroism represents the phenomenon of different 
amount of field absorbtion for the left and right circularly polarized waves as they pass 
through the chiral media. The constitutive relations for chiral metamaterial [8] are 



given below, 



D = eE + «H 



(1) 



B = ^/H - ikE 



(2) 



The behavior of propagation and radiation of the field in a chiral medium has been 
investigated by several authors [8-10]. 

Left-handed materials are those materials in which electric, magnetic and the wave 
vector follow the left hand rule. Vaselago introduced the concept of negative refraction 
for the materials obeying left hand rule [11]. In chiral medium, negative real parts of 
the permittivity and permeability lead an isotropic chiral medium to exhibit circular 
dichroism that is reverse with respect to that exhibited by an identical medium but 
with positive real parts of permittivity and permeability. These left handed materi- 
als exhibit some interesting properties such as backward wave and double negative 
parameters [12]. Negative refraction can also be achieved by using the chiral meta- 
materials [13-16]. In chiral matematerials strong enough chirality produces negative 
refractive index for one of the eigenwaves [17-18]. 

In electromagnetics, Lakhtakia introduced the concept of nihility metamaterial 
for such a material whose permittivity and permeability approaches to zero [19]. Con- 
stitutive relations for nihility metamaterial are, 



Lakhtakia showed that propagation is not possible in nihility metamaterial. Later, 
Tretyakov et al. [20] extended this concept of nihility for the isotropic chiral medium. 



D = 



(3) 



B = 



(4) 



A special case of chiral metamaterial is called chiral nihility metamaterial for which at 
certain value of frequency, real parts of permittivity and permeability both are simul- 
taneously zero. That is, e — > 0, fi — > and k 7^ at nihility frequency. Constitutive 
relations for chiral nihility metamaterials become as, 

D = kH (5) 
B = -i«E (6) 

Chiral nihility metamaterial has two wavenumbers of equal magnitude and opposite 
signs. When a dielectric-chiral nihility half space is excited by an oblique incident 
linearly polarized plane wave, two circularly polarized plane waves having opposite 
handedness are produced. One of them is forward wave while the other is backward 
wave which yields phenomenon of negative refraction in chiral nihility metamaterial. 

If a perfect electric conductor (PEC) interface placed in chiral nihility metamate- 
rial is excited by a forward plane wave. Only one reflected wave parallel to the incident 
wave is produced, i.e., negative reflection happens. The reflected wave is backward 
wave which cancels the incident forward wave to produce zero power propagation. 
Interestingly, if a chiral nihility slab backed by PEC interface is excited by linearly 
polarized plane wave, effect of each eigenwave in chiral nihility after reflection from 
PEC interface is canceled leading to a situation as if front face of the geometry is 
PEC and there is no chiral slab [21]. A PEC waveguide coated with chiral nihility 
metamaterial confines propagation of power within un-coated region of the guide [22] . 

In electromagnetics, Brewster angle is defined as angle of incidence for which re- 
flection power is zero. In all introductory text books on electromagnetic theory [23,24], 



reflection of vertically polarized (parallel polarized or transverse magnetic (TM)) plane 
wave and horizontally polarized (perpendicular polarized or transverse electric (TE)) 
plane wave from a dielectric-dielectric half space boundary is usually presented to ex- 
plain the concept of Brewster angle. It is mentioned that vertically polarized wave 
experiences zero reflection at one incident angle whereas for all angles, reflection is 
non-zero for horizontal polarization. This means that if dielectric-dielectric half space 
is excited by a wave, having both the vertical and horizontal polarizations, at the 
Brewster angle the reflected wave will be linearly polarized with horizontal polariza- 
tion only. The angle of incidence that allows total reflection of power from a planar 
dielectric-dielectric interface is known as critical angle. Critical angle exists only for 
perpendicular polarization if the wave propagates from a denser dielectric medium to 
a rare dielectric. 

In the same context, reflection of plane wave from an interface of an achiral dielec- 
tric and a chiral/chiral nihility metamaterial was discussed by Qiu et al. [25]. Their 
study reveals that for certain values of constitutive parameters, the results similar and 
opposite to that of conventional dielectric-dielectric interface can also be achieved, 
that is, existence of zero reflection for only parallel polarization in the dielectric-chiral 
case. They also showed that, for certain values it is possible to have no-Brewster angle 
for both polarizations and total reflection can be achieved for a wide range of incident 
angles. Power corresponding to the electromagnetic waves associated with a planar in- 
terface of two chiral and/or chiral nihility metamaterials was studied by Faiz et al. [26] 
when it is excited by a plane wave. Using numerical results, interesting characteristics 
such as complete power transmission/rejeection and band pass/band reject filter were 
observed. 



1.1. GO Field and Maslov's Method 

Geometrical optics (GO) is powerful tool for the study of wave motion at high 
frequencies [27-29], however it fails at the the caustics. A caustic is a region where 
area of a ray tube is zero and hence introduces a mathematical singularity when ever 
intensity per unit area is needed to be counted though physically this is not the case. 
Physically the field is always finite at caustics and caustics are of practical interest 
in many applications including defense and medical sciences. Maslov proposed an 
alternative method to find the fields in the caustic region [30]. Maslov's method 
combines the simplicity of asymptotic ray theory and the generality of the Fourier 
transform method. This is achieved by representing the GO fields in terms of mixed 
coordinates consisting of wave vector coordinates and space coordinates. Maslov's 
method are applied to study the field near the caustics of the focusing systems by 
many authors [31-40]. Focusing from a PEC cylindrical and spherical reflectors using 
Maslov's method was first treated by Hongo and Ji [33]. Faryad and Naqvi extended 
this work for reflector coated with chiral metamaterial [39]. Illahi and Naqvi studied 
the focusing from chiral nihility coated PEC and perfect electromagnetic conductor 
(PEMC) cylindrical reflector [41]. In present work, one interest is also to unify all 
these cases so that these become special cases of the geometry, under investigation. 

Consider three dimensional vector wave equation, in cartesian coordinates r = 
(x, y, z), describing the field in medium having wavenumber ho 

V 2 U(r) + k%U(r) = 

Expressing the solution U(r) of wave equation in terms of well known Luneberg-Kline 
series yields Eikonal equation for phase s(r) and transport equation for amplitude. For 



homogeneous and lossless medium, Eikonal equation reduces to Hamiltonian equation 

as 



where p is the wave vector. The solution of Hamiltonian equation is given below 

X = £+p x t, V = V+Pyt, Z = C+p z t, Px=Px () , Py=Py , Pz = Pz 

where (£, 77, () and (p Xo , p yo , p Zo ) are the initial values of cartesian coordinates (x, y, z) 
and wave vector coordinates (p x ,p y ,p z ) respectively and t is parameter along the ray. 
Finite field around caustic may be obtained using following expression obtained from 
Maslov's method [34,40] 



so is the initial phase. Above equation provides uniform solution around the caustic 
region. 

In this paper, transmitted and reflected powers from two parallel layers of chiral 
and/or chiral nihility metamaterials with respect to angle of incidence and chiralities 
of the layers are analyzed. Limiting cases of permittivity, with value of chirality equal 
to zero, are also taken into account and their interpretation are given. Focusing from a 
large size paraboloidal reflector, composed of two layers of chiral and/or chiral nihility 
metamaterials, is studied by utilizing the Maslov's method to give the remedy of GO 
which fails at caustic. 




p = Vs 




exp(-ifcoOo + t - x(p x ,p y , z)p x - y(p x ,p y , z)p y + xp x + yp y ))dp x dp. 



■y 



2.0. Formulation 

Consider a geometry consisting of two infinite parallel layers of chiral metama- 
terials placed in air as shown in Fig. 2.1. Physical width of the two layers is taken 
as di and efe, respectively. Constitutive parameters for the two layers are denoted by 
(ei, /zi, «i) and (e2, ^2,^2), respectively. Air in which layers are placed have constitu- 
tive parameters eo and hq. 

The above mentioned geometry is excited by an oblique incident linearly polar- 
ized electromagnetic plane wave. $i, 9 r and 9 t are the angle of incidence, reflection 
and transmission with respect to the walls of layers, respectively. As chiral medium 
supports two circularly polarized eigen waves having different wavenumbers so fields 
in each layer are written as linear combination of left circularly polarized (LCP) and 
right circularly polarized (RCP) plane waves. Both LCP and RCP waves propagate in 
forward as well as backward direction in each layer. The wavenumbers of both eigen 
waves in a chiral media are different. Refractive indices, wavenumbers and intrinsic 
impedance of the medium for first chiral layer are, 

n l(R,L) = V^l^rl ± «1 (7) 

h(R,L) = w (y/eifii ± «i) (8) 




where n\R and km are refractive index and wavenumber for RCP wave, respectively. 
Quantity e r \ = ei/eo is the relative permittivity of the medium occupying first layer. 



medium 2 




Figure 2.1. Fields representation inside and outside the layers. 



For second chiral layer, these quantities are, 



n 2(R,L) = V /e ^2^r2 ± ^2 
h(R,L) = U (V e 2/U2 ± K 2 ) 



nil 




(10) 

(ii) 

(12) 



Wavenumber and intrinsic impedance of the host medium are given below, 



k =u y/^ e 



(13) 
(14) 



Expressions for the incident and reflected fields and the field transmitted after 
passing through both the layers are, 

E inc = [E i u(a z sm9 i + a x cos 4 ) + E i ±a v ] exp[—iko(z cos $i -xsin/9;)] (15) 



H inc = l/r) {-E i± (a z sm6i + a x cos6i) + E^a y ] exp{-ik (z cos 6 t -xsin^)] (16) 
E re f = [E r \\ (— a z sin9 r + a x cos 9 r ) + E r ±a y ] exp[iko(z cos 9 r + x sin6> r )] (17) 
H re f = l/i]Q[E r ±(—a z sin9 r + a x cos 9 r ) — E r ^a y ] exp[iko(z cos 6 r + x sin6* r )] (18) 
E tra = [E t \\(a. z sm9 t + a x cos9 t ) + E t ±a y ] exp[-ik (z cos 9 t -xsin9 t )] (19) 
H tm = -l/r]o[E t _ L (a z sm9 t + a x cos9 t ) - E t \\a y \ exp[-ik (z cos 9 t -xsm9 t )} (20) 

In above equations, E r \\, E r ±, E t \\, and E t ± are unknown coefficients corresponding 
to parallel and perpendicular components of the respective fields. 

Expressions of fields inside each chiral layer can be written as, 

E p+ = E fL [(a x cos9 pL + a z sm9 pL ) + ia y ] exp[-ik pL (z cos 9 pL -xsm9 pL )} 

+ E f R [(a x cos 9 pR + a z sin 9 pR ) - ia y ] exp[-ik pR (z cos9 pR - x sin9 pR )] (21) 

H p+ = -(i/rj p )Ef L [(a x cos 9 pL + a 2 sin 9 pL ) + ia y ] 

x exp[-ik pL (zcos9 pL - xsm9 pL )} + (i/i] p )E fR [(a x cos9 pR + a z sin9 pR ) - ia y ] 
x exp[-ik pR (z cos 9 pR ] - x sm[9 pR )] (22) 

E p _ = E bL [(-a z sm9 pL + a x cos9 pL ) + ia y ] exp[ik pL (z cos 9 pL +a;sin^p L )] 

+ E bR [(-a z sin 9 pR + a x cos 9 pR ) - ia y ] exp[ik pR (z cos 9 pR + x sm9 pR )] (23) 

H p _ = -(i/r] p )E bL [(-a z sm9 pL +a x cos9 pL ) + ia y ] 

x exp [i k pL { z cos 9 pL + x sin 9 pL )\ + (i/r] p )E bR [(-a z sin9 pR + a x cos9 pR ) - ia y ] 
x exp[ik pR (zcos9 pR + xsm9 pR )] (24) 

where k p L and k pR represent the wavenumbers for LCP and RCP waves, respectively. 
Value of p, as 1 and 2, describes fields for first layer and the second layer, respectively. 



Subscripts b and / are used to represent forward and backward waves, respectively. 
In above equations, -E/l, EfR, E^l, and E^r are unknown coefficients. Tangential 
components of fields are continuous at interfaces located at z = 0, d±, di+d 2 and these 
boundary conditions are used to find the unknown coefficients. Snell's law is used to 
determine the relation among the angles of incidence, reflection and transmission by 
using following relations, 

9 r = 9i 

/c sm/9; = k 1{ R, L ) sm6» 1(RjL) = fc 2 (fl,L) ^^2(r,l) = k t sm6 t 

where ko and k t denote the wavenumbers for medium before and after the layers and are 
equal when these two layers are placed in air. Using transfer matrix method (TMM) 
discussed in [42], three matching matrices each relating the fields at an interface are 
obtained. Two 4x4 propagation matrices, linking fields between two consecutive in- 
terfaces are obtained. Numerical results for powers corresponding to parallel and 
perpendicular components of transmitted and reflected fields are plotted. In section 2, 
it is assumed that optical width of each layer is Ao/4, where Ao is wave length in air 
corresponding to the operating frequency. 

2.1. Results and Discussion 

In this section, numerical results for powers corresponding to parallel and perpen- 
dicular components of the reflected and transmitted fields are presented and analyzed. 
Four different cases, summarized below, are considered for this purpose. 

i. Both layers are of chiral metamaterials (c-c) 



ii. First layer is of chiral nihility and second is of chiral metamaterials (cn-c) 

iii. First layer is of chiral whereas second is of chiral nihility metamaterials (c-cn) 

iv. Both layers are of chiral nihility metamaterials (cn-cn). 

Results corresponding to above cases are shown by Figure 2.2 to Figure 2.9. 

In addition to the above, the limiting cases of permittivity of second layer keeping 
k 2 = are also considered for the discussion. These cases are, 

a. First layer is of chiral metamaterial whereas permittivity of second layer is taken 
very large with chirality k 2 equal to zero (c-PEC) 

b. First layer is of chiral nihility metamaterial whereas permittivity of second layer 
is taken very large with chirality k 2 equal to zero (cn-PEC) 

c. First layer is of chiral metamaterial whereas permittivity of second layer is taken 
very small with chirality k 2 equal to zero (c-n) 

d. First layer is of chiral nihility metamaterial whereas permittivity of second layer 
is taken very small with chirality k 2 equal to zero (cn-n) 

Results corresponding to above situations are shown by Figure 2.10 to Figure 2.13. | 
and _L signs are used to represent power for parallel and perpendicular components, 
respectively. For example: c — cy and c — cj_ represent parallel and perpendicular 
components of power for chiral-chiral cases, respectively. 

In Figure 2.2, powers corresponding to parallel and perpendicular components of 
reflected fields versus angle of incidence are shown. It is assumed that the intrinsic 
impedances of metamaterials filling both layers are same. Figure contains plots for four 



different cases, i.e., c-c, cn-c, c-cn, and cn-cn. Brewster angle is defined as the angle of 
incidence for which reflection of power, both parallel and perpendicular components, 
is zero. Reflected powers for c-c case are zero for 6>j < 23° and this gives a range 
of brewster angles. In the case of cn-cn, a range of Brewster angles exists for 9i < 
17.5°. After this range, power for parallel component increases whereas power for 
perpendicular component remains negligible for the entire range of incident angles. In 
the case of cn-c, parallel component is nonzero everywhere and near grazing incidence 
almost total reflection of power happens whereas power for perpendicular component 
remains negligible through the entire range of incident angles. It may be noted that 
range of Brewster angles exist only for c-c and cn-cn cases and range of angles yielding 
almost total reflection exists only for parallel component of cn-c case. For c-cn case, 
power only for perpendicular component is zero for 9i < 27.5° and corresponding 
parallel component is nonzero everywhere. 

Figure 2.3 describes behavior of the parallel and perpendicular components of 
transmitted powers for impedance matching of both layers. For c-c case, power for 
parallel and perpendicular components have same initial amplitude, i.e., 0.5. Consid- 
ering cn-c case, power of both parallel and perpendicular components of transmitted 
power follow each other. In c-cn case, perpendicular component is dominant over par- 
allel component whereas for cn-cn case, it is observed that power of parallel component 
is dominant over the perpendicular component. 



ii and x reflected power vs angle of incidence 
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Figure 2.2. Behavior of the reflected power for parallel and perpendicular field com- 
ponents. Ki = .25, K2 = -75, T]i=T]2 = 2 

Figure 2.4 shows the behavior of reflected power for layers having different intrinsic 
impedances. For c-c and cn-c cases, no Brewster angle exists and power for parallel 
component is dominant over corresponding perpendicular component. For c-c case, 
reflected power only for perpendicular component is zero for 9i < 34.3° and 9{ > 85°. 
For cn-c case, when 9i > 22°, power for parallel component of reflected field starts 
increasing and after 9i = 80°, almost total reflection is observed because perpendicular 
component has insignificant reflected power. Perpendicular component has very low 
value between #^=34.3° and 0j=78° and is zero elsewhere. Both c-cn and cn-cn cases 
have same range of Brewster angles, i.e., 9i < 11.46°. 

Figure 2.5 is about transmitted power in case of impedance mismatch. For c-c 
case, both components of power have equal contribution at 9i = 31.5°. For cn-c and 
cn-cn cases, power for parallel component has higher values and follows the power for 
corresponding perpendicular component. 





nents. K\= K2 = .25, 771 = 1 and 772 = 2 




ii and J- reflected power vs chirality/C[ 




Figure 2.6. Behavior of reflected power for parallel and perpendicular field compo- 
nents. K2 = .25, 9i = 7r/4 and 771 =772 = 2 



Figure 2.6 shows plots of power for impedance matching case versus chirality of 
the first layer. For the purpose of analysis chirality values are considered between 
and 5. For cn-c case, almost total power reflection in terms of parallel component 
occurs at very low values of chirality and perpendicular component is almost zero for 
considered range of chirality. For c-cn case, zero reflection of parallel power compo- 
nent is obatined at chirality value k\ = 1.9 whereas at K\ = 1.8,2.2,4 perpendicular 
component of reflected power is zero. In cn-cn case, oscillating behavior is observed 
for both components with periodic behavior for cn-cn^. At K\ = 0, 2,4 perpendicular 
component of power becomes zero and at K\ = 1,3, 5 parallel component of power is 
zero. At K\ = 0, total reflection occurs. Zeros of parallel and perpendicular compo- 
nents of power at specific values of chirality can be used to get or avoid particular 
component of polarization. 

Figure 2.7 gives transmission behavior of power. It is noted that both components 
for each case contain either zeros or values approaching to zero. These findings may 
be used get or avoid particular component of polarization. In Figure 2.8, no Brewster 
angle exist for all four cases except for cn-c case where both components are negligibly 
small at ki = 1.9. It is also noted that perpendicular components of power for cn-cn 
and cn-c are negligibly small. Conclusions similar to Figure 2.7 may be drawn from 
Figure 2.9. 
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Figure 2.7. Behavior of transmitted power for parallel and perpendicular field com- 
ponents. K2 = -25, 9i = 7r/4, T]i=T]2 = 2 

Figure 2.10 shows the plots for reflected power assuming very high value of per- 
mittivity of second layer. This situation may be considered as a chiral layer backed by 
perfect electric conductor interface. In general reflected power has both components, 
i.e., both c — PEC|| and c — PEC^ are nonzero. But for certain values of chirality 
of first layer, reflected field contains only co-polarized component. Contribution of 
parallel component is very large as compared to the perpendicular component of re- 
flected power. For cn-PEC case, complete reflection in terms of parallel component is 
observed. Figure 2.11 shows zero transmission of power for c — PEC and cn — PEC 
cases. 

Another interesting case is shown in Figure 2.12 where second layer is assumed 
of nihility metamaterial. This has been achieved by taking «2 = and setting very 
low value of permittivity of second layer. This situation may be seen as a chiral layer 



C-C 




backed by nihility interface. For c-n case, almost total reflection in terms of parallel 
component of the reflected power is observed. For the case of cn-n, perpendicular 
component of reflected power has also significant contribution for certain range of 
chirality. It is noted that for certain range of chirality, overlap of cn-n and c-n cases 
happen. Figure 2.13 shows plots of transmitted power for c-n and cn-n cases and 
no transmission of power is obtained. It may be noted that conclusions drawn from 
Figure 2.10 to Figure 2.13 agree with published work. 




ii and x transmitted power vs chirality^ 2 
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Figure 2.9. Behavior of transmitted power for parallel and perpendicular field com- 
ponents. Ki = .25, 771 =772 = 2, and 6i = tv/4 
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Figure 2.10. Behavior of reflected power for parallel and perpendicular field compo- 
nents. 9i = 7r/4, 7/i = 2, e 2 = 10 5 , k 2 = 
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Figure 2.11. Behavior of transmitted power for parallel and perpendicular field com- 
ponents. 6i = 7r/4, ?7i = 2, e 2 = 10 5 , k 2 = 
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Figure 2.12. Behavior of reflected power for parallel and perpendicular field compo- 
nents. 6i = 7r/4, rji = 2, e 2 = 10~ 5 , k 2 = 
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Figure 2.13. Behavior of transmitted power for parallel and perpendicular field com- 
ponents. 6i = ir/4, rji = 2, e 2 = 10~ 5 , k 2 = 



Focusing of electromagnetic waves from a paraboloidal reflector composed of chi- 
ral and/or chiral nihility metamaterial using the Maslov's method is studied in next 
section. PEC and nihility backed chiral/chiral nihility paraboloidal reflector are also 
discussed. 

3.0 Evaluation of Finite Field Around Focus 

Considered a bilayer paraboloidal reflector placed in air as shown in Figure 3.0. 
The equation for the surface of a paraboloidal reflector is given by, 

( = m,v) = f-^ = f-^l (26) 

where (£, rj, () are the Cartesian coordinate of the point on the surface of paraboloidal 
reflector. / is the focal length of the paraboloidal reflector and p 2 = £ 2 + rj 2 . 

The incident field traveling along z axis is expressed as, 

Ej = 3l x exp (-ik z) (27) 

The incident plane wave make an angle a with surface normal, where surface normal 
is given as, 

a n = sin a cos 7a x + sin a sin ja y + cos aa z (28) 
where a and 7 are given as, 

cosa=^====, (30) 
tan 7 = I (31) 




Figure 3.0. Reflection from a paraboloidal reflector placed in air 
The wave reflected from paraboloidal reflector is given below, 

E r = E(r ) exp[iko(x sin 2a cos 7 + y sin 2a sin 7 + z cos 2a)] (32) 

Rectangular components of field E(r ), at the surface of reflector, are given below, 

E xQ = B± sin 2 7 — Bn cos 2 7 cos 2a (33) 

E y0 = — cos7sin7(S|| cos 2a + B±) (34) 

E z o = B» sin 2a cos 7 (35) 



where B\\ and B± are parallel and perpendicular components of GO reflected field 



from planar layers of chiral and/or chiral nihility layers. The field expression in polar 
coordinates and valid around the focal point is given below [40], 



x E r0 tanoiexp[— ik (2f — r sin6>sin2aicos(</> — 7) — r cos 6 cos 2a)]dad r y (36) 
The upper limit of integration with respect to a is taken as, 



where D is the height of the paraboloidal reflector from horizontal axis. 

3.1 Numerical Results and Discussion 

The field behavior around the focal region of a paraboloidal reflector are obtained 
by solving the equation, obtained using Maslov's method. Numerical results are pre- 
sented with the variation of different parameters, i.e., chirality, relative permittivity 
and thickness of the layer. Chiral nihility is introduced as limiting case of chiral 
medium where value of relative permittivity and relative permeability of the chiral 
medium is assumed as 10~ 5 , keeping nonzero chirality whereas for nihility material, 
the value of chirality parameter is taken equal to zero. In all the simulations, param- 
eters (k ,f,H) are taken as (1,100, 7r/4) . For all cases, permeability of each chiral 
layer is taken as unity. All plots deal with behavior of |U(r)| versus kz. 

3.1.1 Chiral-Chiral Paraboloidal Reflector 

From Figure 3.2 and Figure 3.3, it is noted that by increasing the permittivity of 
each layer of the paraboloidal reflector, increases the strength of electric field around 





the focus. Comparison of Figures 3.2 and 3.3, also shows that strength of field around 
focus is more dependent on the permittivity of second layer. 

In Figure 3.4 and Figure 3.5 , effects of variation of thickness of each layer of 
paraboloidal reflector are presented. It is noted that increase of thickness of first or 
second layer also increases the strength of field around focus. It is noticed that first 
layer has more impact on the field strength as compared to the variation of thickness 
of the second layer. 

Figures 3.6 and 3.7 show the effect of change of chirality parameter of each layer 
while keeping all the other parameters fix. It is observed that increase in chirality of 
each layer also strengthen the field around focal point of the paraboloidal reflector. 
Change in chirality of first layer produces greater change in reflected field strength 
around the focal point as compared to chirality variation for second layer. 

For the case of c-c paraboloidal reflector, some interesting features have been 
studied that are permittivity, chirality and thickness of each chiral layer is directly 
proportional to the reflector's field strength at focal point. Moreover it has been 
observed that, first and second layers show different effects against the variation of a 
particular parameter while keeping all the other parameters fix. 

3.1.2 Chiral-PEC and Chiral-Nihility Paraboloidal Reflector 

In this section two special cases for paraboloidal reflector are discussed one by one: 
paraboloidal reflector composed of c-PEC and c-n cases. Figure 3.8 shows reflected 
field strength for three different values of permittivity of the first layer when second 
layer is considered as PEC (very large value of permittivity, i.e., 10 5 ). It is noted that 



increase in permittivity of first layer also increases the field strength at the focal point. 
Figure 3.9 shows that increase in the chirality of first layer correspondingly increases 
the field strength at the focus. In Figures 3.10 and 3.11, second layer is considered as 
nihility (by assuming very small values of permittivity and permeability, i.e., 10 -5 ). 
From Figures 3.10 and 3.11, it is observed that by increasing the values of permittivity 
and chirality, the field strength around the focal point also increases correspondingly. 

4.0 Conclusions 

Characteristics of reflected and transmitted powers from a structure of two par- 
allel layers filled with chiral and/or chiral- nihility metamaterials are studied. Special 
cases of PEC and nihility backed chiral/chiral nihility layer are also discussed. It is 
observed that the behavior of reflected or transmitted powers strongly depend on the 
angle of incidence and values of constitutive parameters describing the matamaterials. 
For some specific values and ranges of incident angles, total reflection of power has 
been observed hence structure can yield range of Brewster angles. Parallel or per- 
pendicular component of reflected/transmitted power may be selected. That is, either 
co-polarized or cross polarized field is produced. Similarly for some values of chirality, 
total reflection/ transmission of power is also observed. This arrangement of parallel 
layers can also be used as a power divider. 

Arrangement of two parallel layers mentioned above is utilized to study the behav- 
ior of field around the focal region of a large size paraboloidal reflector. Paraboloidal 
reflector, composed of c-c, c-PEC and c-n layers is discussed for this purpose. PEC 
and nihility backed cases give higher values of reflected field around the focus as com- 
pared to composed of chiral-chiral layers. It is concluded that increase in the relative 



permittivity of each layer also increases field strength around the focus. An increase 
in the thickness and chirality parameter of the layer also have the same effect, i.e., 
increase in values of these parameters yields correspondingly higher focusing of the 
field. It is also noted that different parameters, mentioned above, affect the strength 
of the field differently. Moreover in the case of c-c paraboloidal reflector, variation of 
a parameter for first and second layer yields different reflected field strength when all 
other parameters are kept constant. In this case, change in permittivity of second layer 
produces more variation in the reflected field strength around the focal point as com- 
pared to the first layer, whereas for change in chirality and thickness, field strength of 
the reflector is more sensitive for the first layer. Hence from above discussion strength 
of the reflected field around the paraboloidal reflector can be controlled by studying 
the effect of different factors. Study of focusing field strength is of great importance 
in defence and other electronic systems. This study provides an investigation to de- 
sign an efficient reflector coated with materials having different material properties. 
This whole discussion suggests to get control over the focusing of the reflected field by 
bringing variations in the permittivity, thickness or the chirality of the corresponding 
layer. 
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Figure 3.2. Field behavior around focus for different values of permittivity of first 
layer: d\ = di = c-c case. 
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Figure 3.3. Field behavior around focus for different values of permittivity of second 



layer: d\ = d,2 = 4 , c-c case. 





Figure 3.5. Field behavior around focus for different values of thickness of second 
layer: d\ = ^f, c-c case. 
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Figure 3.6. Field behavior around focus for different values of chirality of first layer: 



d± = d,2 = -T-, c-c case. 



\U(r)\ 

40 r 



30 



ki = .25 
e rl = e r2 = 1 



20 



k 2 = .90 
k 2 = .75 
k 2 = .50 




-20 



-10 







10 



20 



kz 



Figure 3.7. Field behavior around focus for different values of chirality parameter of 
second layer: d\ = (1% = c-c case. 




Figure 3.9. Field behavior around focus for different values of chirality of first layer: 
d x = d 2 = ^f, e r2 = 1 x 10 5 , c-PEC case. 




Figure 3.11. Field behavior around focus for different values of permittivity of first 
layer: d\ = ^f, e r 2 = 1 x 10~ 5 , k 2 = 0, c-n case. 



References 



1. D. F. Arago, "Sur une modification remarquable qu' eprouvent les rayons lu- 
mineux dans leur passage a travers certains corps diaphanes, et sur quelques 
autres nouveaux phenomnnes d'optique", Mem. Inst. 1, 93-134 (1811). 

2. J. B. Biot, "Ph6nomenes de polarisation successive, observs dans des fluides ho- 
mogenes", Bull. Soc. Philomath., 190-192 (1815). 

3. A. Fresnel, "Memoire sur la double refraction que les rayons lumineux eprouvent 
en traversant les aiguilles de cristal de roche suivant des directions parall61es a 
1'axe", Oeuvres 1, 731- 751 (1822). 

4. J. B. Biot, "M6moire sur la polarisation circulaire et sur ses applications a la 
chimie organique", Mem. Acad. Sci. 13, 39-175 (1835). 

5. K. F. Lindman, "Ober eine durch ein isotropes system von spiralformigen res- 
onatoren erzeugte rotationspolarisation der elektromagnetischen wellen", Ann. 
Phys. 63, 621-644 (1920). 

6. K. F. Lindman, "Uber die durch- ein aktives raumgitter erzeugte rotationspolar- 
isation der elektromagnetischen wellen", Ann. Phys. 69, 270-284 (1922). 

7. D. L. Jaggard, A. R. Mickelson, C. H. Papas, "On electromagnetic waves in chiral 
media", Appl. Phys. 18, 211-216 (1979). 

8. I. V. Lindell, A. H. Sihvola, S. A. Tretyakov, A. J. Viitanen, "Electromagnetic 
waves in chiral and bi-isotropic media", Artech House, Boston (1994). 

9. S. Bassiri, C. H. Papas, and N. Engheta , "Electromagnetic wave propagation 



through a dielectric-chiral interface and through a chiral slab", J. Opt. Soc. of 
Am. A, 5, 1450-1459 (1988). 

10. A. Lakhtakia, "Beltrami fields in chiral media", World Scientific, Singapore(1994). 

11. V. G. Vaselago, "The electrodynamics of substances with simultaneously negative 
values of permittivity and permeability", Sov. Phys. Usp. 10, 509-514 (1968). 

12. C. W. Qiu, H. Y. Yao, S. Zouhdi, L. W. Li, and M. S. Leong, "On the constitutive 
relations of G-chiral media and the possibility to realize negative-index media", 
Microwave Opt. Technol. Lett. 48, 2534-2538 (2006). 

13. S. Tretyakov, A. Sihvola, and L. Jylh, "Backward- wave regime and negative re- 
fraction in chiral composites", Photonics Nanostruct. Fundam. Appl. 3, 107-115 
(2005). 

14. J. B. Pendry, A. J. Holden, W. J. Stewart, and I. Youngs, "Extremely low fre- 
quency plasmons in metallic mesostructures" Phys. Rev. Lett. 76, 4773-4776 
(1996). 

15. J. B. Pendry, "A chiral route to negative refraction" Science 306, 1353-1355 (2004). 

16. C. W. Qiu, H. Y. Yao, L. W. Li, S. Zouhdi, and T. S. Yeo, "Routes to left-handed 
materials by magnetoelectric couplings" Phys. Rev. B 75, 245214 (2007). 

17. S. Zhang, Y. S. Park, J. Li, X. Lu, W. Zhang, and X. Zhang, "Negative refractive 
index in chiral metamaterials" , Phys. Rev. Lett. 102, 023901 (2009). 

18. T. G. Mackay, A. Lakhtakia, Simultaneous negative- and positive-phase-velocity 
propagation in an isotropic chiral medium, Microwave and Optical Technology 
Letters 49, 1245-1246 (2007). 



19. A. Lakhtakia, "An electromagnetic trinity from negative permittivity and negative 
permeability", Int. J. Inf. and Mil. Wav., 22, 1731-1734 (2001) 

20. S. Tretyakov, I. Ncfedov, A. H. Sihvola, S. Maslovki, C. Simovski," Waves and 
energy in chiral nihility", Journal of Electromagnetic Waves and Applications 17, 
695-706 (2003). 

21. Q. A. Naqvi, "Planar slab of chiral nihility metamaterial backed by fractional 
dual/pemc interface", Progress In Electromagnetics Research PIER 85, 381-391 
(2008). 

22. M. A. Baqir, A. A. Syed and Q. A. Naqvi, "Electromagnetic fields in a circular 
waveguide containing chiral nihility metamaterial" Progress In Electromagnetics 
Research M 16, 85-93 (2011). 

23. C. A. Balanis, "Advanced engineering electromagnetics", John Willey and Sons, 
2nd Edition (2012). 

24. D. K. Cheng, "Fields and wave electromagnetics", Addison- Wesley, Newyork 
(1989). 

25. C. W. Qiu, N. Burokur, S. Zouhdi, and L. W. Li, "Chiral nihility effects on energy 
flow in chiral materials", J. Opt. Soc. of Am., 25, 55-63 (2008). 

26. F. Ahmad, S. N. Ali, A. A. Syed, and Q. A. Naqvi, "Chiral and/or chiral nihility 
interfaces: parametric dependence, power tunneling and rejection" Progress in 
Electromagnetics Research M 23, 167-180 (2012). 

27. L. B. Felson, Hybrid formulation of wave Propagation and scattering, Nato ASI 
Series, Martinus Nijho, Dordrecht, The Netherlands, (1984). 



28. G. A. Dechamps, "Ray techniques in electromagnetics," Proc. IEEE, 60, 1022- 
1035 (1972). 

29. C. H. Chapman and R. Drummond, "Body wave seismograms in inhomogeneous 
media using Maslov asymptotic theory," Bull. Seismol., Soc. Am., 72, 277-317, 
(1982). 

30. V. P. Maslov, "Perturbation theory and asymptotic method" , Gos. Moskov., 
Univ., Moscow, (1965) (in Russian). (Translated into Japanese by Ouchi et al., 
Iwanami, Tokyo, 1976). 

31. A. Ghaffar, Q. A. Naqvi, and K. Hongo, "Analysis of the fields in three dimensional 
Cassegrain system", Progress In Electromagnetics Research PIER 72, 215-240, 
(2007). 

32. Y. Ji, and K. Hongo, "Analysis of electromagnetic waves refracted by a spherical 
dielectric interface by Maslov's method", J. Opt. Soc. Am. A, 8, 541-548 ( 1991). 

33. Y. Ji, and K. Hongo, "Field in the focal region of a dielectric spherical by Maslov's 
method", J. Opt. Soc. Am. A, 8, 1721-1728 (1991). 

34. K. Hongo, , Y. Ji, and E. Nakajima, "High frequency expression for the field in 
the caustic region of a reflector using Maslov's method", Radio Sci. 21, 911-919 
(1986). 

35. K. Hongo, and Y. Ji, "High frequency expression for the field in the caustic region 
of a cylindrical reflector using Maslov's method", Radio Sci. 22, 357-366 (1987). 

36. K. Hongo, and Y. Ji, "Study of the field around the focal region of spherical 
reflector antenna by Maslovs method", IEEE Trans. Antennas Propagat., 36, 



592-598 (1988). 

37. R. W. Ziolkowski, and G. A. Deschamps, "Asymptotic evaluation of high fre- 
quency field near a caustic: an introduction to Maslov's method" , Radio Sci. 19, 
1001-1025 (1984). 

38. M. Faryad, and Q. A. Naqvi, "High frequency expression for the field in the caustic 
region of cylindrical reflector placed in chiral medium" , Progress In Electromag- 
netics Research PIER 76, 153-182 (2007). 

39. M. Faryad, and Q. A. Naqvi, "High frequency expression for the field in the caustic 
region of a parabolic reflector coated with isotropic chiral medium," Journal of 
Electromagnetic Waves and Applications 22, 965-986 (2008). 

40. T. Rahim, , M. J. Mughal, Q. A. Naqvi, and M. Faryad, "Focal region field of a 
paraboloidal reflector coated with isotropic chiral medium" Progress In Electro- 
magnetics Research PIER 94, 351-366 (2009). 

41. A. Illahi and Q. A. Naqvi "Study of focusing of electromagnetic waves reflected 
by a PEMC backed chiral nihility reflector using Maslov's method", Journal of 
Electromagnetic Waves and Applications 23, 863-873 (2009). 

42. C. Sabah, S. Uckun "Mirrors with chiral slabs" Journal of Optoelectronics and 
Advanced Materials 8, 1918-1924 (2006). 



